

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
lim_angl=80;
id=1;
noise=.03;
image='phantom';
maxit=10;
initialGuessStr='zero';
constThreshValue=10e-8;
imgSize=256;
mode='';%'smooth_lin_';
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


angles=(-lim_angl/2):(lim_angl/2);
truncate=true;



switch id
    case 1
        identifier=['limitedAngle_outerFading_' num2str(lim_angl)];
    case 2
        identifier=['limitedAngle_innerFading_' num2str(lim_angl)];
    case 3
        identifier=['limitedAngle_fixedFading_3_' num2str(2/3*lim_angl) '_' num2str(lim_angl)];
end

if(~strcmp(mode,''))
    identifier={{identifier},{[0 1]},{[mode num2str(2/3*lim_angl) '_' num2str(lim_angl)]}};      
end

switch image
    case 'phantom'
        img=phantom(imgSize);
    case 'lena'
        img=double(imread('lena512.pgm'))/255;
        sd=512/imgSize;
        if(sd>1)
            img=sample_down(img,sd);
        end
end

Y=radon(img,angles);

if(noise>0)
    Y=addNoise(Y,noise);
end


out1=my_iterative_soft_thresholding(Y,angles,imgSize,maxit,initialGuessStr,constThreshValue,identifier)
%out2=iterative_soft_thresholding(Y,angles,truncate,imgSize,maxit,initialGuessStr,constThreshValue)